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Abstract 



A new method of constructing a weak coupling expansion of two dimensional 
. models with an unbroken continuous symmetry is developed. The method 

is based on an analogy with the abelian XY model, respects the Mermin- Wagner 
(MW) theorem and uses a link representation of the partition and correlation 
functions. An expansion of the free energy and of the correlation functions at 
small temperatures is performed and first order coefficients are calculated explic- 
itly. They are shown to coincide with the results of the conventional perturbation 
theory. We discuss an applicability of our method to analysis of uniformity of 
the low-temperature expansion. 



1 Introduction 

Since two-dimensional models with a continuous global symmetry were recognized to 
be asymptotically free [TJ they became a famous laboratory for testing many ideas 
and methods before applying them to more complicated gauge theories. In this paper 
we follow this common way and present a method, different from the conventional 
perturbation theory (PT), which allows to investigate these models in the weak coupling 
region. The conventional PT is one of the main technical tools of the modern physics. 
In spite of the belief of the most of community that this method gives the correct 
asymptotic expansion of such theories as 4D QCD or 2D spin systems with continuous 
global symmetry, the recent discussion of this problem has shown that it is rather 

far from unambiguous solution. Indeed, for the PT to be applicable it is necessary 
that a system under consideration would possess a well ordered ground state. In 
two dimensional models, like the 0(iV)-sigma models the MW theorem guarantees 
the absence of such a state in the thermodynamic limit (TL) however small coupling 
constant is |7|. Then, the usual argument in support of the PT is that locally the system 
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is ordered and the PT is not supposed to be used for the calculation of long-distance 
observables. On the other hand it should reproduce the correct behaviour of the fixed- 
distance correlations as well as all thermodynamical functions which can be expressed 
via short-range observables. The example of ID models shows that even this is not 
always true ||, so why should one believe in correctness of the conventional PT in 2D! 
In fact, the only way to justify the PT is to prove that it gives the correct asymptotic 
expansion of nonperturbatively defined models in the TL. Now, it was shown in || that 
PT results in 2D nonabelian models depend on the boundary conditions (BC) used to 
reach the TL. This result could potentially imply that the low-temperature limit and 
the TL do not commute in nonabelian models. Actually, the main argumentation of 
|| regarding the failure of the PT expansion is based on the fact that the conventional 
PT is an expansion around a broken vacuum, i.e. the state which simply does not 
exist in the TL of 2D models. According to ||, the ground state of these systems can 
be described through special configurations - the so-called gas of superinstantons (SI) 
and the correct expansion should take into account these saddle points. At the present 
stage it is rather unclear how one could construct an expansion in the SI background. 
Fortunately, there exists other, more eligible way to construct the low-temperature 
expansion which respects the MW theorem and is apriori not the expansion around 
the broken vacuum. We develop this method on an example of the 2D SU(N) x SU(N) 
principal chiral model whose partition function (PF) is given by 

^ReTvU x U+ +n , (1) 

x,n 

where U x G SU(N), DU X is the invariant measure and we impose the periodic BC. 
The basic idea is the following. As was rigorously proven, the conventional PT gives 
an asymptotic expansion which is uniform in the volume for the abelian XY model |§ . 
One of the basic theorems which underlies the proof states that the following inequality 
holds in the 3D XY model 

< exp( v /^A(0 1 .)) > < C , (2) 

where C is /3-independent and A(<p + 2tt) = <fi. <p x is an angle parametrizing the action 
of the XY model S = J2xn cos (0a; — <fix+n)- It follows that at large f3 the Gibbs measure 
is concentrated around <j) x ~ providing a possibility to construct an expansion around 
<j) x = 0. This inequality is not true in 2D in the thermodynamic limit because of the 
MW theorem, however the authors of M prove the same inequality for the link angle, 
i.e. 

< exp(y^A(0i)) > < C , (f)i = <p x - <p x+n , (3) 

where the expectation value refers to infinite volume limit. Thus, in 2D the Gibbs 
measure at large j3 is concentrated around (pi m and the asymptotic series can be 
constructed expanding the action in powers of (pi. In the abelian case such an expansion 
is equivalent to the expansion around <f) x = because i) the action depends only on 
the difference <p x — <f) x + n and ii) the integration measure is flat, DU X = d(f) x . 
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In 2D nonabelian models, again because of the MW theorem, one has to expect in 
the TL something alike to @, namely 



< exp( v //5argA(Tr^)) > < C , V l = U x U+ +n . (4) 

Despite there is not a rigorous proof of (|]), that such (or similar) inequality holds in 2D 
nonabelian models is intuitively clear and should follow from the chessboard estimates 
TUfl and from the Dobrushin-Shlosman proof of the MW theorem M which shows that 



spin configurations are distributed uniformly in the group space in the TL. Namely, 
the probability that Tr(VJ — I) < — £ is bounded by 

p(£) < 0(e- b ^) , -> oo , (5) 

if the volume is sufficiently large, b is a constant. Thus, until £ < 0((^/J3)~ 1 ) is not 
satisfied, all configurations are exponentially suppressed. This is equivalent to the 
statement that the Gibbs measure at large (3 is concentrated around VJ ~ I, therefore 
(|) or its analog holds. In what follows it is assumed that (ffl) is correct, hence V\ = I 
is the only saddle point for the invariant integrals^. Thus, the correct asymptotic 
expansion, if exists, should be given via an expansion around Vi — I, similarly to the 
abelian model. If the conventional PT gives the correct asymptotics, it must reproduce 
the series obtained expanding around Vi = I. However, neither i) nor ii) holds in the 
nonabelian models, therefore it is far from obvious that two expansions indeed coincide. 
Let us parametrize V/ = exp(iu;;) and U x = exp(iu x ). Consider the following expansion 

The standard PT states that to calculate the asymptotic expansion one has to re- 
expand this series as 

wj = u x - uj x+n + -^kj^^ » ( 7 ) 

where uu^ are to be calculated from the definition U x U x+n = exp(iuj[). This is pre- 
sumably true in a finite volume where one can fix appropriate BC (like the Dirichlet 
ones), or to break down the global symmetry by fixing the global gauge (on the periodic 
lattice). Then, making f3 sufficiently large one forces all the spin matrices to fluctuate 
around U x ~ /, therefore the substitution ([?[) is justified. We do not see how this 
procedure could be justified when the volume increases and fluctuations of U x spread 
up over the whole group space. In other words, it is not clear why in the series 



ui — uj x — uj x+n + Lo 

k=l 



(*) 



4 It follows already from (||) . What is important in is a factor \/]5, otherwise the very possibility 
of the expansion in 1//3 becomes problematic. 
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the term u\ is suppressed as /5~ 1//2 relatively to the term u\ . It is only (|6|) which 
remains correct in the large volume limit and takes into account all the fluctuations 
contributing at a given order of the low-temperature expansion. 

It is a purpose of the present paper to develop an expansion around V\ — I aiming 
to calculate the asymptotic series for nonabelian models. First of all, one has to give 
a precise mathematical meaning to the expansion (P). It is done in the next section. 



2 Link representation for the partition and corre- 
lation functions 

To construct an expansion of the Gibbs measure and the correlation functions using 
(§) we use the so-called link representation for the partition and correlation functions. 
First, we make a change of variables V\ = U x U^ +n in ([I]). PF becomes 



Y[ dVi exp 



/3^Re TrVi + \nJ(V) 



where the Jacobian J(V) is given by |T 



J(V)= JU dU -Ii EdrXrty+U^ 



n 



l€p 



(9) 



flp is a product over all plaquettes of 2D lattice, the sum over r is sum over all 
representations of SU(N), d r = Xr(I) is the dimension of r-th representation. The 
SU(N) character Xr depends on a product of the link matrices Vi along a closed path 
(plaquette in our case): 



1] V = V n (x)V m (x + n)V+(x + m)V+{x) . 



(10) 



lep 



The expression J2r d r Xr(J\ie P V) is the SU(N) delta-function which reflects the fact 
that the product of U x U+ +n around plaquette equals / (original model has L 2 degrees 
of freedom, L 2 is a number of sites; since a number of links on the 2D periodic lattice 
is 2L 2 , the Jacobian must generate L 2 constraints)^. The solution of the constraint 



Uv = i 

lep 

is a pure gauge Vi = U x U x+n , so that two forms of the PF are exactly equivalent. 
The corresponding representation for the abelian XY model reads 



\\ d(f>i exp 



P cos <frl 



I 



114 



p ■ 



(12) 



j p 



5 Strictly speaking, on the periodic lattice one has to constraint two holonomy operators, i.e. closed 
paths winding around the whole lattice. We do not expect such global constraints to influence the TL 
in 2D (see Discussion). 
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where the Jacobian is given by the periodic delta-function 

oo 

J p = E elT(t>P ' = + ^rn{x + n) - (f) n (x + m) - 4> m (x + n) . (13) 

r=— oo 

Consider two-point correlation function 

F{x,y) =< Tr U x U+> , (14) 

where the expectation value refers to the ensemble defined in ([I]). Let C xy be some 
path connecting points x and y. Inserting the unity U Z U£ in every site z G C xy one 
gets 

T(x,y) =< Tr [J {U m U+J >=< Tr J] W t > , (15) 

where W\ = V\ if along the path C xy the link / goes in the positive direction and 
Wi = V^, otherwise. The expectation value in fllST) refers now to the ensemble defined 
in (|). Obviously, it does not depend on the path C xy which can be deformed for 
example to the shortest path between sites x and y. 

In this representation the series (^j) acquires a well defined meaning, therefore the 
expansion of the action, of the invariant measure, etc. can be done. 



3 XY model: Weak coupling expansion of the free 
energy 

In this section we prove that for the abelian XY model the large-/? expansion in the link 
representation gives the same results as the conventional PT in the thermodynamic 
limit. We consider only the free energy but the generalization for the correlation 
functions is straightforward. 

The first step is a standard one, i.e. we rescale — > -4= and make an expansion 



exp 



(3 cos 



VP 



exp 



where l\ + 2l 2 + ... + klj~ = k and 



ai 1 



fc=i 



h,-,h 



\fc+l 



(,2(fc+l) 

(2fc + 2)! ' 



..ai 



h\...l k \ 



(16) 



(17) 



In addition to this perturbation one has to extend the integration region to infinity. 
We do not treat this second perturbation, as usually supposing that all the corrections 
from this perturbation go down exponentially with (3 (in the abelian case it can be 
proven rigorously |J). It is more convenient now to go to a dual lattice identifying 
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plaquettes of the original lattice with its center, i.e. p — *■ x. Let I = (x; n) be a link on 
the dual lattice. Introducing sources hi for the link field, one then finds 



Z XY {f3 » 1) = e 



/3DL 2 -L 2 \n/3 



n 



00 1 . ( d 2 

k 



dhl 



M(hi) . 



Coefficients Ak are defined in fll6|) and fllTl). The generating functional M(hi) is given 
by 



E 



n« 

1 



n exp 



-E0? + *E^= 
2 i 1 VP 



r x - r. 



+ E <M* 



(19) 



Sums over representations are treated by the Poisson resummation formula. The Gaus- 
sian ensemble appears in terms of the fluctuations of r-fields. The integral over zero 
mode of r-field is not Gaussian and leads to a delta-function in the Poisson formula. 
Thus, the zero mode decouples from the expansion. All the corrections to the integrals 
over representations in the Poisson formula fall down exponentially with j3 so that the 
generating functional becomes 



M(hi 



exp 



1 



E hiGiyhi 



u> 



(20) 



where we have introduced the following function which we term "link" Green function 



Gw — 2oi t i' — G X)X i — G x+ri)X > +n i + G XyX / +n i + G x+n ^ 
G XyX i is a "standard" Green function on the periodic lattice 



I L-l e ^L kn (x-x') n 
Gx,x' — Y2 E 



L 2 k % D - cos fk n 



(21) 



(22) 



Normalization is such that Gu = 1. As far as we could check, the expression 
reproduces the well known asymptotics of the free energy of the XY model in two 
dimensions. For example, the first order coefficient of the free energy being expressed 
in Gw reads 



64L 2 



E^ 2 



1 

32 



(23) 



Let us add some comments. In the standard expansion to avoid the zero mode problem 
one has to fix appropriate BC, like Dirichlet ones or to fix a global gauge if one works 



on the periodic lattice |12j| . In the present scheme the zero mode decouples automati- 
cally due to using £7(1) delta-function which takes into account the periodicity of the 
integrand in link angles. More important observation is that it is allowed to take the 
TL already in the formula (|T8D: since the generating functional depends only on the 
link Green function which is infrared finite, the uniformity of the expansion in the 
volume follows immediately. This is a direct consequence of the fact that the Gibbs 
measure of the XY model is a function of the gradient (pi only. 
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4 Weak coupling expansion in the 577(2) model 



We turn now to the nonabelian models. As the simplest example we analyze the SU(2) 
principal chiral model. The method developed here has straightforward generalization 
to arbitrary SU(N) or SO(N) model and we shall present it elsewhere. 



4.1 Representation for the partition function 

We take the standard form for the £77(2) link matrix which is the most suitable for 
the weak coupling expansion 

V t = exp[z<7 fc w fc (Z)] , (24) 
where a k , k = 1,2,3 are Pauli matrices. Let us define 



and similarly 



w n 



5>*(0 



L k 



1/2 



1/2 



where u>f.{p) is a plaquette angle defined as 



lep 



exp 



(25) 



(26) 



(27) 



The exact relation between link angles 10^(1) and the plaquette angle ouk(p) is given 
in the Appendix B. Then, the partition function @ can be exactly rewritten to the 
following form appropriate for the weak coupling expansion 



n 

i 



sin 2 Wi 



Y[du k (l) 



cxp 



2(3 cos Wi 



n 



w x 



n e 

x m(x)=— 



rfafc(x) exp 

fc 



sin W 7 ^ 



(28) 



where we have introduced auxiliary field a^{x) and 

.1/2 



a(x) 



L k 



(29) 



The representation for the Jacobian (]71|) used here can be interpreted as an analog of 
the U (1) periodic delta-function (|I3|): this is the SU (2) delta-function which is periodic 
with respect to a length of the vector uj(x). We give a derivation of the partition 
function (p5[) in the Appendix A. In rewriting the final formula of that derivation ([73]) , 
we went over to the dual lattice identifying the links of the original lattice with dual 
links and the original plaquettes with dual sites located in the center of the original 
plaquettes. 
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4.2 General expansion 



To perform the weak coupling expansion we proceed in a standard way, i.e. first we 
make the substitution 

u k (l) - (2[3)- l / 2 u k (l) , a k (x) - (2f3) 1/2 Mx) (30) 

and then expand the integrand of ( p8|) in powers of fluctuations of the link fields. We 
would like to give here some technical details of the expansion which could be useful 
for a future use. It is straightforward to get the following power series: 

1. Action 



exp 



2(3 cos 



Wi 



exp 



2/5 --W) 2 



af...a£ 



1 + EK0-' E ilL . y 



(31) 



fc=i 



where Zi + 2/2 + ••• + /c/fc = k and 



Ofc = (-1) 



fc+i v n 



W, 



2(fe+l) 



(2k + 2)\ 



(32) 



2. Invariant measure 



sin 2 Wi 00 (— l) fc fc 1 

~WT = 1 + g W ^ ' ^ = 5 (2n + l)!(2ife - 2n + 1)! • (33) 

3. Contribution from Jacobian I 



TT/ OO T r)2fc — 1 

^ = 1 + Et^^, ^ = 2^— |5 2fc |, (34) 
sin 14^ ^ (2(3 ) k (2k)\ 



where B 2k are Bernoulli numbers. 
4. Contribution from Jacobian II 



a k (x)u k (x) = a k (x) 

exp -iJ2a k (x)u k (x 

00 

i + E 



4 0) (*) + E 



^ n) (*) 



i (2/9)"/ 2 



(35) 



exp 



9=1 



HI! 



00 ,.,(")/ 

E^W E 

fc n= 



- (2/3)W2 
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Using the relations ([T5|)-(|80|) one can calculate oo^ k n \x) up to an arbitrary order in n. 
In particular, uf\x) is given by (see Fig.l for our notations of dual links) 



uJ k \x) = UJ k (l 3 ) + uj k (l A ) - uj k (h) - UJ k (l 2 ) . 



We shall use an obvious property 



^a k (x)uf\x) = ^2u) k (l)[a k (x + n) -a k (x)} , I = (x;n) . 



(36) 



(37) 



Introducing now the external sources h k (l) coupled to the link field oj k {l) and s k (x) 
coupled to the auxiliary field a k (x) and adjusting the definitions 



u k (l) 







a k (x) 







dh k (l) ' ' ' ds k (x) 
we get finally the following formal weak coupling expansion for the PF (|28|) 

%c k wr 



Z = C(2P)Z(0,0)l[ 



n 

X 

where 



Ju 



1 (2/3) 



2k 



l + E(2/3)- fc E 

k=l h,--,h 
(-<)« 



Zl Ih. 

gj — a k 
hl..l k \ 



1 (2/3)' 



1+ S * 



C(/3) = exp 



00 ^ n) (z) 



? a * (ar) S(2^ 



2/ 3L 2 -^L 2 ln/3 



M(/i, s) 



(38) 



(39) 



(40) 



As usually, one has to put h k = s k = after taking all the derivatives. M(h, s) is a 
generating functional which we study in the next subsection. It is obvious that the 
ground state satisfies 

<(J k °\x)Y>=0, p= 1,2,... , (41) 
precisely like the abelian model. 



4.3 Generating functional and zero modes 

Here we are going to study the generating functional M(h, s) given by 

Z(h,s) 



and 
Z(h, s) 



Y[ da 



k(x 



x,k 
oo 



M(h,s) 



Y[du k (l) exp 

l,k 



Z(0,0) 



(42) 



-uj 2 k (l) - iu k (l)[a k (x + n) - a k (x)] 



E ex p 

m(x)=— oo 



2Tciyf2f3j2 m ( x ) a ( x ) + E^(0M0 +E Q; fc( :r ) s fc( :r ) 

x Lk x,k 



(43) 
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As in the abelian case we expect that integrals over zero modes of the auxiliary field 
are not Gaussian and should lead to some constraint on the sums over m x . To see this, 
we put h k = s k = and integrate out the link fields. Partition function becomes 



Z(0,0) 



E 

m(x)=— oo 



k{x)exp 



x,k 



—a k (x)G x x ,a k (x') + 2mJ2f3m(x)a(x 



(44) 



with G X:X i given in ( P2|) and sum over repeating indices is understood here and in what 
follows. It is clear that in this case the zero modes should be controlled via integration 
over the radial component of the vector a(x). To see this, we change to the spherical 
coordinates and treat only a constant mode in the angle variables (since the zero mode 
problem could only arise from this configuration). One has 



00 roo 

Z(0,0)~ / Y[a 2 x da x exp 



m x =— 00 
00 



a x G x x ,a x i + 2mJ2j3m x a 



r /m r > 



+ 0(mt 



(45) 



and we used the notation a x for the radial component of the vector a(x). Since the 
zero mode of the radial component in the x space is 



a{p = 0) 



'^E(E«^)) 2 ) 



1/2 



one has to omit the zero mode from the Green function in each term of the sum over 
k in the integrand of (f4"4"p. As in the abelian case the integration over zero modes 
produces delta-function in (|45"D. Since however only m x = for all x contribute to the 
asymptotics of the free energy and fixed distance correlation function (other values of 
m x being exponentially suppressed) we have to put m x = omitting at the same time 
all zero modes. Calculating resulting Gaussian integrals we come to 



M(h, s) = exp 



1 i 1 

-s k (x)G X}X >s k (x') + -s k (x)Di(x)h k {l) + -h k (l)G u ,h k (l') 



where Gu> was introduced in ([H]) and 

Di(x) = G x>x > - G x+riiX/ , l = (x,n) . 
From fl4l)| ) one can deduce the following simple rules 

< u k (l)uj n (l f ) >= -^-Giv , < a k {x)a n {x') >= ^G^ , 

-i < u k {l)a n (x ) >= ^yDi(x') . 



(46) 



(47) 



(48) 



We describe some simple properties of the functions Gw and D[(x') in the Appendix 
C. The expansion fl40|), representation ( [46D for the generating functional and rules ( ^8|) 
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are main formulas of this section which allow to calculate the weak coupling expansion 
of both the free energy and any short- distance observable. Let us now comment on the 
infrared fmitness of the expansion. It follows from the representation for the generating 
functional that all expectation values of the link fields are expressed only via the 
link Green functions Gw and Di(x) which are infrared finite by construction. All 
combinations of auxiliary fields which contain odd overall powers of the fields are 
expressed only via Di(x) and, therefore are also infrared finite. However, even powers 
include G XjX > and the infrared finitness is not provided automatically. In particular, it 
means that unlike XY model we are not allowed to take the TL at this stage. 

Our last comment concerns the partition function (JPj). We believe it can be re- 
garded as an analog of the corresponding expression in the XY model, i.e. this is 
a nonabelian analog of the so-called "spin-wave-vortex" representation for the par- 
tition function. One can see that the nonabelian model is not factorized into three 
abelian components and is periodic in the length of the vector a(x) rather than in its 
components ctk(x). 

4.4 First order coefficient of the correlation function 

As the simplest example we would like to calculate the first order coefficient of the 
correlation function flT5|). Expanding ([15]) in 1/(3 one has 

r(*,y) = i-^<£(l> fc (o) > +o(/r 2 ) = 1 - E g ui + o([3- 2 ), (49) 

where is a path dual to the path C xy , i.e consisting of the dual links which are 
orthogonal to the original links I, I' G C xy . The form of Gw ensures independence of 
r(x, y) of a choice of the path C xy . After some algebra it is easy to get the result 

!■(*, y) = 1 - ±D( X - y) , D (x) = ^ g ^'f^ • « * », (50) 
which coincides with the result of the conventional PT. 



5 First order coefficient of the free energy 

The main result of our study is the first order coefficient of the SU(2) free energy 

F = ^ln Z = 2(3 -|ln (3 + -i^C 1 + 0(/T 2 ) • (51) 
There are four contributions at this order to C 1 

C 1 = Cac + Cmeas + @J1 + Cj2 • (52) 
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Contribution from the action fl3T|) is given by 



1 1 5 2 5 

C ac = H G U = ■ ( 53 ) 



2L 2 ac 128L 2 Y 64 



Contribution from the measure fl33| ) is given by 

^Las = -^E G « = -^ ( 54 ) 

Contribution from third brackets in fl4"0|) is proportional to [^i°' > ] 2 and equals 0, because 
of (PI). There are two contributions from the expansion of the Jacobian (0). The 
first one is given by the expectation value of the operator — i Y, x J2l=i a k(x)oo^\x). 
uj^\x) is given in the Appendix B (|83D . From the form of the generating functional 
([46|) one can see that the expectation value of this operator depends only on link Green 
functions Gw and Di(x'). One gets after long but straightforward algebra^ 

1 ^4eeaw 



2L 2 " 4L2 , , 



1 4 



/ i=i 

$ih(Gi 3 i 4 + 2Gi s i 2 + 2Gm 2 ) + 5u 2 (G[ 3 i 4 — G^ — G^iJ + 
SihiGi^ + Gi 4 i 2 — Gi ± i 2 ) — du^Gi^ + 2Gi 3 i 2 + Gi ± i 2 ) ) , (55) 

where links U are defined in Appendix C (see Fig.l), I = (x,n). In terms of standard 
/^-functions defined in Appendix C the result reads 

1 iC), = ~[6 - 2D{2, 0) - D(l, 1)] = i + i- . (56) 



2L 2 Ji 4 l v > y 2 2tt 

The second term is given by the operator < | J2l=i a k(x)uJ^\x)j >. u^\x) is 
given in the Appendix B (|32"D. In terms of Green functions it reads 

\^ 2 = -i(Q (1) + Q (2) ), (57) 



2L 2 Jl 16 

4 4 

Q {1) = ^J2J2J2 G xA G k,vG h v }l - G w > ; G/.^ ; ) , (5* 

a;, 2' i<j i'<j' 



Q {2) ^EEE^^^WA.M + ^^^WAix') 



Gi^i' Dp {x)Di.{x') - G lhV D v (x)D li {x')) . (59) 



3 Our previous version suffered from incorrect sign in this term which led to wrong final result. 
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Link li (I'j,) refers to one of four links attached to a given site x (x'). We performed both 
analytical and numerical studies of last expressions. Details are given in the Appendix 
C. We find 

QW = \ (60) 

and 

Q (2) = 1 + - . (61) 

7T 

Collecting all coefficients we finally obtain 

' C 1 = £ , (62) 



2L 2 64 

which coincides with the result of the conventional PT 



6 Discussion 

In this paper we propose to use an invariant link formulation to investigate some prop- 
erties of 2D models in the weak coupling region. We argued that this approach is more 
suitable for calculation of asymptotic expansions of invariant functions in cases when 
the Mermin- Wagner theorem forbids spontaneous symmetry breaking in the thermo- 
dynamic limit. We have found that both in the abelian XY model and in non-abelian 
SU(2) model our results for the first order coefficients of the free energy and correlation 
function agree with the standard PT expansion in the TL. In our next paper [[13] we 
show that the second order coefficient of the correlation function also agrees with the 
conventional PT. 

Now we can return to the question raised in the Introduction, namely whether 
conventional PT gives uniform asymptotic expansion for non-abelian models. It is well 



known for a long time that it is no so in one- dimensional non-abelian models. In |13| 
we address this question in the link formulation and show that non-uniformity in this 
case originates from the expansion of holonomy operator which imposes certain global 
condition on the configurations of link matrices. While this global condition itself 
vanishes in true TL, it survives large volume limit if the low-temperature expansion 
is performed in a finite volume. In 2D there are two such operators which we did 
not consider in the present paper. The reason being that the contributions from these 
holonomy operators are suppressed as 0(1/ L). Since in 2D one could encounter only 
logarithmic divergences it is rather unlikely that holonomies are relevant for the TL. In 
this sense there is no similarity between one and two dimensional models. Nevertheless, 
one cannot exclude the possibility of non-uniformity of the low-temperature expansion 
arising from the remainder to the PT series ||. This problem is extremely hard to 
resolve by means of the standard approaches, see for instance |J. Contrary, in the 
link formulation we are able to calculate the exponential remainder at a given order of 
the low-temperature expansion. Therefore, the problem of the infrared finitness of the 
remainder can be addressed explicitly. Such investigations are presently in progress. 



13 



Acknowledgements. 

We are grateful to J. Polonyi who found time to go through many details of the 
expansion presented here and for many encouraging discussions. We would like to 
thank V. Miransky and V. Gusynin for interesting discussions and healthy critics on 
different stages of this work. Our special thanks to B. Rusakov for the explanation of 



his paper \'\ 1] regarding the calculation of the Jacobian J(V). 



References 

[1] E. Brezin, J. Zinn- Justin, Phys.Rey. B14 (1976) 3110; E. Brezin, J. Zinn- Justin, 
J.C. LeGuillou, Phys.Rey. D14 (1976) 2615; A.M. Polyakov, Phys.Lett. B59 (1975) 
79. 

[2] A. Patrascioiu, E. Seiler, Phys.Rev.Lett. 74 (1995) 1920. 
[3] A. Patrascioiu, E. Seiler, Phys.Rev.Lett. 74 (1995) 1924. 
[4] F. David, Phys.Rev.Lett. 75 (1995) 2626. 

[5] F. Niedermayer, M. Niedermaier and P. Weisz, Questionable and Unquestionable 
in the Perturbation Theory of non-Abelian Models, [hep-lat/9612002| . 

[6] A. Patrascioiu, E. Seiler, Phys.Rev. D57 (1998) 1394. 

[7] N.D. Mermin, H. Wagner, Phys.Rev.Lett. 22 (1966) 1133; R. Dobrushin, S. Shlos- 
man, Com.Math.Phys. 42 (1975) 31. 

[8] Y. Brihaye, P. Rossi, Nucl.Phys. B235 (1984) 226. 

[9] J. Bricmont, J.-R. Fontaine, J.L. Lebowitz, E.H. Lieb, and T. Spencer, 
Comm.Math.Phys. 78 (1981) 545. 

[10] J. Frohlich, R. Israel, E.H. Lieb, B. Simon, Comm.Math.Phys. 62 (1978) 1. 

[11] B. Rusakov, Phys.Lett. B398 (1997) 331. 

[12] P. Hasenfratz, Phys.Lett. B141 (1984) 385. 

[13] O. Borisenko, V. Kushnir, in preparation. 



14 



7 Appendix 

7.1 A: Representation for SU(2) partition function 

We start from the following partition function in the link representation 



Z = J n^exp 



P E TiV i 



n w 



Using the formula for SU(2) characters 

Xr(4>) 



sin(2r + 1) 



sm i 



with r = 0, 1/2, 1, 3/2, we can write down the Jacobian J(V) © in the form 



J P (v)= E 



sm r, 



r v =-oo Sm( Pp 



(63) 



(64) 



(65) 



where <p p is some plaquette angle. Now r p takes only integer values. Let us parametrize 
the SU (2) link matrices as 



Vi = exp ia k uo k (l) = cos W\ + ia k uo k (l) 



sin Wi 
Wi 



(66) 



where we used W\ defined in (|2"5D. One gets in this parametrization the following 
expressions: 



1. Action 



2. Invariant measure 



3. Jacobian 



-TrV, = cos Wi 



sin 2 Wi 



J P (Y)= E 



Y[duj k (l) . 



sin r p W p 
smW n 



(67) 
(68) 

(69) 



where W p is defined in (^6]). Substituting (|6TD- (j69|) into (|63| ) we get 
"sin 2 Wi 



n 



exp 



n e 



sin r^W 7 , 



P r p =-oo P sinWp 



(70) 



To get an expression for the Jacobian convenient for the large-/? expansion we use the 
equality 



~ smrW W 
r hj sinW 7 ~ siniy„ ; _ x ., , 



oo . 

E / ]Jda k exp 



-i E a k^k + 2iTima 

k 



(71) 
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where we have introduced a = [J2k a k\ • ^° P rove (0) we wr ^ e 



/oo 
\{da k e-^ / 
-OO i J — i 



k k 



ia k t k 



(72) 



sinri 



Integrals over t k are calculated in the spherical coordinates. Taking then F(t 2 
and using the Poisson resummation formula we come to (|7T|) . Substituting ( f7T|) into 
the partition function flTOj ) we finally get 



n 



sin 2 Wi 



k 



exp 



2(3 c °s Wi 



n H * 



p sin W p 



oo - 



P m(p)=— oo 



-i a k (p)uj k (p) + 2%im(p)a(p) 



(73) 



7.2 B: Relation between link and plaquette angles 

Let us introduce the following notations for links of a given plaquette p 

l\ = (x\ n) , l 2 = (x + n; m) , Z 3 = (x + m; n) , U— (x\ m) . (74) 

Then, from the definition (p7|) and ( |T0l) one gets the following formulas relating link 
and plaquette angles 



V p = cos W p + itrukip) — p 



where 



cos VTp = cos Mi cos M 2 + ^ i/f 



fc fe sin Mi sin M 2 



'A. 



M X M 2 



sin W n 



,,,,, w / p k . . sinMi k sinM 2 fc „ „ a sinMisinM 2 
u; fc (p) p = v\ cosM 2 — i/£ cos Mi—- h C vq v\v\ — — 



Mi 



M, 



Mj and i/,- are given by 

cos Mi = cos W(Zi) cos W(Z 2 ) - ^2 ^k{h)^k{h) 

k 

cosM 2 = cos W(h) cosH / (/ 4 ) -^^(^3)^(^4) 



MiM 2 

sinV^(/i)sinV^(/ 2 ) 

W(/i)W(/ 2 ) ' 
sinW(Z 3 )sinW/(Z 4 ) 



(75) 

(76) 
(77) 

(78) 



fc sinMi 
z/f— — = u k (h) cos W(Z 2 ) 
ik/i 



W(l 3 )W(h) 
sinW(h) . ......siii^lj) 



^ + u k (l 2 ) cos W(h 



W(hj 



k smM 2 ,. , TXrn . smW(l 3 ) . TTT/1 

Vo . ; = UJ k (l 3 ) COS W(l 4 ) Tlrn \ +U k {k) C0SW(1 



W{1 2 ) 
smW{lx)sinW{h) 
W{h)W{l 2 ) 



(79) 



Mo 



e kpq u p (h)uj q (h) 



W{h) 
sinW / (/ 3 )siniy(/ 4 ) 
W(h)W(U) 
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Now it is straightforward to calculate the following expansion {jp — > x) 

u k {x) = uf\x) + J k l \x) +4\x) + ... . (80) 
On a dual lattice (see Fig.l) the first coefficients can be written down as 

uJ ( k\ x ) = u k {h) + u k (U) - u k (li) - uj k (l 2 ) , (81) 

J k 1 \x) = -e k ™j:u p (h)cu q (l J ), (82) 

i<j 

J k 2 \x) = \e™ e k ^\x)jz^{k)^ m {k) + 

-(e e +e e )[u r (l 3 )LO m (li)uj n (l2) + Ur{h)Um{h)u>n{h) - 

U r {h)^m{k)^n{k) - UJ r (h)^m(h)^n(k)] • (83) 

7.3 C: Coefficients and 

We define D-function as 

D(xux 2 ) = — £ 9 v2 —WZ 0, (84) 



and Gn as 



l l 



Go = T2 T. 2 T 2 CQS i« k , fc^O. (85) 

^ fc„=0 Z ^n=l COS T K « 

The normalization is chosen such that Gu = 1. All .D-functions can be expressed in 
terms of D(l, 0) and D(l, 1). In the TL one has 

D(l,0) = i 13(1,1) = -, 

I 7T 

D(2, 0) = 2(2£>(1, 0) - D(l, 1)) = 2 - - , D(2, 1) = 2D(1, 1) - D(l, 0) = - - ~ , 

7T 7T Z 

D(3, 0) = 4(13(2, 0) - D(l, 1)) + 0) = i + 8(1 - -) . (86) 

Z 71 

Gq is known to diverge logarithmically in two dimensions and behaves numerically as 
G = a + ai InL + o(l) , a » 0.0974006 , a x w 0.318331 « - . (87) 

7T 

We also need some representation for the function Gw in the momentum space. One 
finds from the definition (|21f) 

F(n = n') = 2(1 — cos —k p ) ,n p , 

±j 

F(n ^ n') = (1 - e^ fc ")(l - e"^') , (88) 
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x I, 



Figure 1: Plaquette of original lattice and links of dual lattice as they enter sums for 
and Q^ 2 \ Link is determined by point x and a positive direction, e.g. Z3 = (x — ni). 



2 27T. 



where we have denoted 

f(k) = 2-Y,cos^k n . (89) 

n=l ^ 

Using this representation it is easy to prove the following "orthogonality" relations for 
the functions G\y and Di(x) which are very useful for calculation of lattice sums 

Y,G lb G hV = 2G lv , (90) 

b 

"£D b (x)G w = 0, (91) 

b 

Y,D h {x)D h {x') = 2G x>x , . (92) 

b 

To visualize the summation over indices in (0) and (|59| ) we depicted dual links 
in Fig.l. 

7.3.1 Q<V 

We divide into two pieces 

Q« = B X + B 2 , (93) 

where B\ includes first and second powers of Green functions G XyX i whereas B 2 consists 
of terms with only third powers of G x>x '. We thus have for B\ from (p8|) 

1 44 



2L 2 

x,x' i<j v <y 



Performing all summations and using Eq. (|8T)|) to express all D functions appearing in 
the last equation in terms of D(l, 0) and D(l, 1) we get 

B x = 4D(1, 1) - 4D(1, 0) - 4G Q . (95) 

For we adduce an expression in the momentum space 

2 2U ^ + P) ' 1 J 
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where 



B(k 1 n , kl) = (zi + z 2 + z{ + z 2 + z 2 z\ + z 1 z 2 )(z 4 + z 3 + z 3 + z% + z 4 z 3 + z 3 zl) + (97) 
2(«i + z*)(^4 + z* 4 ) - 2z 1 z*z 4 z 3 - 2ziz 2 (z 3 z 4 )* - (ziz 2 + {ziZ2y){zzZ A + (z 3 z 4 )*) - 

4(4? - zHzif - 244(4 + 4) - ^Mi4f + (4) 2 ) + (4 - - (4) 2 ) 

and we have introduced 

1199 Ait 
z x = 1 - e'* , z 2 = 1 - e»* , z 3 = 1 - e l * , z 4 = 1 - e ip * ; p* n = —k\ . (98) 



The result of summation can be expressed in terms of D functions given in (SB) 



B 2 = 5£>(1,0)- 4£>(1, 1) + 4G . (99) 
One sees that all divergences exactly cancel. The final result is 

= B 1 + B 2 = D(1,0) = ^. (100) 

We performed numerical check of our result for B 2 . The function was calculated for 
lattice sizes L e [10 — 120] and fitted to the form 

B 2 = a + a 1 \\iL + a 2 -^ + -^ . (101) 

From the result = ~ and from Eq.(|87D one concludes that coefficients a and a 4 
should be equal to 

a = 0.34312331 , a x = 1.273324 . 

Our fit gives 

a = 0.343476 , a x = 1.27327 . 

7.3.2 

As before we divide into two pieces 

Q® = K t + K 2 , (102) 

where K\ includes second powers of Green functions, K 2 is cubic in G xx >. Ki is given 
by 

a;,a;' «<j i'<j' 

5 W , ; ^>)A 3 (^) - ^A^A^')) • (103) 
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Calculations are rather lengthy nevertheless straightforward. We get 

8 



K x =20D(l,0r + 4(D(2,0)-D(l,0))(2D(l,l)-3D(l,0) + D(2,0)) = 8-- . (104) 

71 

K 2 in the momentum space can be written as 

Nil"*™ 

where 

^(^n> fc n) = (*L + ^2 + ^3 + Z 4 ~ Z t Z 3 - Z 2 Z A )[(z 1 + Z 2 f (z* z + zj) 

-(21 + z'l + z 2 + z 2 )(z 3 + z 4 )] 
-2 Zl z* 4 [2(z 2 + l)(z* + z* A - z*zl) + z x {z 2 + z 4 - z 2 z 4 )] . (106) 

Notations are as in Eq.(K8J). All sums can be done analytically resulting in[] 

16 

K 2 = — -7. (107) 

7T 

Finally we obtain 

Q( 2 ) = - + 1 « 3.5464791 . (108) 

7T 



7 This analytical expression agrees well with our previous numerical result = 3.5466309 
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